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Abstract 

This paper completes the argument that lattice simulations of partially quenched 
QCD can provide quantitative information about QCD itself, with the aid of par- 
tially quenched chiral perturbation theory. A barrier to doing this has been the 
inclusion of ^q, the partially quenched generalization of the rj' , in previous calcu- 
lations in the partially quenched effective theory. This invalidates the low energy 
perturbative expansion, gives rise to many new unknown parameters, and makes it 
impossible to reliably calculate the relation between the partially quenched theory 
and low energy QCD. We show that it is straightforward and natural to formulate 
partially quenched chiral perturbation theory without ^q, and that the resulting 
theory contains the effective theory for QCD without the r]' . We also show that 
previous results, obtained including ^q, can be reinterpreted as applying to the 
theory without <I>o- We contrast the situation with that in the quenched effective 
theory, where we explain why it is necessary to include ^q. We also compare the 
derivation of chiral perturbation theory in partially quenched QCD with the stan- 
dard derivation in unquenched QCD. We find that the former cannot be justified as 
rigorously as the latter, because of the absence of a physical Hilbert space. Finally, 
we present an encouraging result: unphysical double poles in certain correlation 
functions in partially quenched chiral perturbation theory can be shown to be a 
property of the underlying theory, given only the symmetries and some plausible 
assumptions. 
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1 Introduction 



Numerical simulations of lattice QCD are hampered by the difficulty of including loops of 
light quarks. This has forced the use of approximations to the fermion determinant: the 
quenched approximation (setting the determinant to a constant), and, more recently, the 
partially quenched (PQ) approximation (including the determinant but with sea quark 
masses different from, and typically larger than, those of the valence quarks). While all 
such simulations correspond to unphysical theories, they are not all equally unphysical. It 
has been argued recently that PQ simulations can be used to obtain physical parameters 
if the quarks are light enough that one can use chiral perturbation theory to describe 
the low energy properties of the theory |^, |^. The only approximation necessary is 
the truncation of chiral perturbation theory. On the other hand, if the sea-quarks are 
too heavy, then partial quenching is an uncontrolled approximation whose results will at 
best be a qualitative guide to those in the physical theory. 

The main purpose of this paper is to complete the theoretical argument justifying 
the use of PQ simulations to obtain physical parameters. The missing ingredient in 
the arguments presented in Refs. |^ concerns the flavor-singlet field, $o; which is the 
generalization of the v( in PQ theories. This field must be kept in the low-energy effective 
theory for quenched QCD, because its correlation functions contain poles at the masses 
of the light pseudo-Goldstone mesons. The same holds true for the PQ theory with 
heavy sea quarks. The need to include $o invahdates the standard chiral power counting 
and introduces additional coupling constants in the chiral Lagrangian. Consequently, 
progress with this theory can only be made by making further assumptions. 

On the other hand, in QCD the i]' is not a pseudo-Goldstone boson, the corresponding 
field need not be included in the chiral Lagrangian, and chiral perturbation theory can 
be developed as a low energy expansion |0, ^, ^. What we demonstrate here is that the 
situation is similar for PQ QCD with sea (and valence) quarks in the chiral regime: It 
can be described by a chiral Lagrangian in which $0 is absent, and for which standard 
power-counting applies. As stressed in Refs. |Q ^, an important corollary is that the 
parameters of the PQ chiral Lagrangian are the same as those in the chiral Lagrangian 
for QCD. This is what is needed to show that PQ simulations can be used to extract 
physical parameters. 

Another result shown here is of a more technical nature. All previous calculations 
using PQ chiral perturbation theory have included the $0 field, and thus suffer from the 
problems described above. Here we show that these problems can be avoided by sending 
the $0 mass parameter, mo, to infinity, because this is mathematically equivalent to 
considering the theory in which $0 is absent. In this way, the results of the previous 
calculations from PQ chiral perturbation theory can be reinterpreted as applying to the 
theory for which the matching to QCD is immediate. This observation also justifies the 
ad- hoc prescriptions for integrating out $0 that were used previously |]^, ||, ^, ^. 

A secondary purpose of this paper is to discuss the theoretical foundations of chiral 
perturbation theory for PQ theories, building on the work of [PH], |TT|. We recall the line of 



reasoning used to construct the chiral Lagrangian for QCD, examine the extent to which it 
applies also to PQ QCD, and point out the gaps that make PQ chiral perturbation theory 
stand on less secure grounds. We are able to show, however, that a signature prediction 
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of PQ chiral effective theories, namely the presence of double pole contributions to flavor 
singlet correlators, can be derived in the underlying theory from the assumption that 
valence flavor non-singlet correlators have single poles. The latter assumption is well 
tested by numerical simulations. 

This paper is organized as follows. In the next section we review the definition of 
PQ theories, and discuss their symmetries. The set of two-point functions containing 
light poles is identified in Section || for both PQ and quenched theories. In Section | 
we construct the PQ effective Lagrangian that reproduces the previously determined 
pole structure, and discuss its theoretical foundations. In Section |^ we explain how 
chiral perturbation theory including $o is equivalent to that excluding $o if one sends 
rriQ oo. We then discuss the pole structure of fiavor-singlet correlators in Section 
and summarize our conclusions in Section |^. 

Five appendices contain technical details. Appendix A describes the true symme- 
tries of the PQ QCD partition function, derives the resulting Ward identities, and shows 
how these are in fact the same as those obtained assuming a "fake" symmetry group. 
Appendix Ij| collects useful results on the diagonal generators of graded Lie groups. \Kp 



pendix gives an alternative argument for why one does not need to include $o in 
PQ chiral perturbation theory. Appendix D| concerns a useful result about the limit 
mo — > oo. Finally, [Appendix E| derives the constraints on neutral pion correlators using 
graded symmetries. 



2 Partially quenched theories and their symmetries 

We consider a theory with sea quarks and A'^v' valence quarks, which is viewed as a 
tool to study an unquenched theory with N quarks. We will shortly consider all quarks, 
sea and valence, to be light, although this is not needed to set up the definition of the 
theory. Clearly, the theories we have in mind are those with N = 2 (treating only the up 
and down quarks as light) and = 3 (also including the strange quark — whose status 
as a light quark is less clear). The number of valence quarks one usually uses is Ny = 2 
(needed to discuss simple meson properties), or Ny = 3 (needed for baryon properties). 

In practice, one obtains PQ results by generating gauge configurations including in 
the weight a determinant representing the effect of sea quarks, and then calculating 
quark propagators on these background gauge fields using masses which are different from 
those of the sea quarks. This can be represented theoretically by Morel's construction 
involving bosonic spin-1/2 ghost fields |jl2|, in which the Euclidean partition function is 

4) (p + m) + (p + rh) tjj 



Z = j Dl'^ijiPipA] exp (^~Sg - j 



, ^ N det (E) + m) , , 

D A exp -gc? J 1 

det (JZ) + m) 

Here 5*^ is the gauge action, ip an + A^-dimensional column vector containing the 
quark fields, with mass matrix m, and tp an A'^v'-dimensional vector of ghost quarks, with 
mass matrix m. For each valence quark there is a ghost quark with the same mass, so 
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that the valence quark determinant cancels precisely against that from the ghost quarks. 
The complete mass matrix is 

^"('o ^^=^^^s{mvi,...,rnvNv','>TT'Si,...,msN]mvi,...,mvNv)- (2) 

In the following, we will take all non-zero entries of M to be real and positive. 

Our discussion will concern the theory obtained in the continuum limit of PQ sim- 
ulations. This allows us to take p as having standard continuum properties: it is an 
anti-hermitian operator which connects left-handed fields to right-handed fields and vice- 
versa. In this way we avoid the issue of how best to discretize fermions, and of possible 
difficulties in simulating odd numbers of dynamical quarks. We simply assume that these 
difficulties have been overcome, and that the lattice simulations are done close enough to 
the continuum limit that eq. (|l|) represents them up to small corrections, suppressed by 
powers of the lattice spacing, which can be extrapolated away.Q Indeed by using overlap 
fermions, or other fermions with an exact chiral symmetry, one can presumably formulate 
the discussion of symmetries at non-zero lattice spacing. 

Correlation functions of quark and ghost fields are defined, as usual, by introducing 
source terms into Z. Note that if one considers correlation functions involving only 
sea-quark fields, one obtains exactly the result of the unquenched theory with N sea 
quarks, because the valence and ghost determinants cancel [|lOl- Furthermore, if one of 



the valence quarks has the same mass as a sea quark, then it can replace that sea quark 
in correlation functions without changing the result . 



To discuss symmetries, it is useful to collect quarks and ghost quarks into a single 
+ 2A^y-dimensional vector, Q, defined by 

Q^ = (^^^^). (3) 
The fermionic part of the action then takes the standard form 

Sf = [ Qip + M)Q (4) 

{QlPQl + QrPQr + QrMQl + QlMQr) , 

where the projections are defined by 

(1±75) ^ ^ 7^ (1 T 75) 
Ql,r = ^ Q , Ql,r = Q — 2 — ■ 

The symmetry group of the action eq. (^, when M —* 0, appears to be the graded 
group 

U{Nv + N\Nv)l®U{Nv + N\Nv)r, (6) 
under which the fields transform in the usual way 

Ql,r{^) — ' Ul,rQl,r{^) ' QlA^) — ' QlA^^^Ir > (7) 



^Some comments on the impact of 0(a) corrections are made in Ref. 
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with Ul,r G U{Nv + A^|A^y). As noted in Ref. |]TT|, however, this is not the correct 
symmetry because the functional integral over the bosonic spin-1/2 fields converges only 
if they are properly constrained (and if the mass matrix is positive definite). As a result of 
the constraint, right- and left-handed fields cannot be rotated independently in the usual 
way. Nevertheless, it turns out that one can proceed as if the symmetry group were (0), 



as long as one only considers small transformations. This is explained in [Appendix A|. In 



particular we show that one obtains the correct vector and axial Ward identities if one 
pretends that the action has the symmetry group eq. (^), rather than the actual symmetry 
group. In what follows we mostly use the "fake" symmetries, since this emphasizes the 
similarities to the development for unquenched QCD, and is usually simpler. We use the 
real symmetries only when global aspects of the symmetry group are important. 

Certain of the transformations in eq. (|^) are anomalous, since they do not leave the 
measure invariant. Removing these requires that we impose the constraint 

sdet Ul = sdet Ur , (8) 

where "sdet" is the invariant "super-determinant" .0 The overall result is that the chiral 
symmetry group of the massless PQ theory can be taken to be the group proposed in 

m 



SU{Nv + N\Nv)l ® SU{Nv + N\Nv)r ® U{l)v • (10) 

Here the U{l)v transformations are common overall phase rotations of the right- and 
left-handed fields.^ 

The form eq. (|T0|) refiects the fact that for the PQ case (i.e. > 0), the U{1) 
factor can be chosen to be a fiavor singlet, and thus to commute with all elements of 
SU{Nv + N\Nv)- In particular, the anomalous symmetry current can be chosen to be 
the fiavor singlet, = Qliil^Q (using the notation of [Appendix B| , eq. ( [115[) ). This has 



the same quantum numbers as the pseudoscalar density (575<5, which is the interpolating 



field associated with $o |T31 — focus of much of the subsequent discussion. 

The symmetry structure is different for the fully quenched theory, A^ = 0. In partic- 
ular, as noted in Ref. [|T3|, the flavor singlet U{1) is part of SU{Ny\Ny), while the U{1) 



with non-unit superdeterminant cannot be chosen to be a flavor singlet, and so does not 
commute with all elements of SU{Ny\Ny). This is explained more fully in [Appendix B[, 



where we collect some results on the generators of graded groups. The net effect is that 
the SU{Ny\Ny) and U{1) factors form, locally, a semi-direct product, and the chiral 
symmetry group can be taken to be 

[^V{Ny\Ny)L ® ^'\]{Ny\Ny)R] K U {l)y . (11) 



^Here the super-dctcrminent is defined by sdet U = expstrlnJJ, with the supertrace being 

str ^ = tvA - tiD , (9) 

the blocks corresponding to the Ny + N quark and Ny ghost coordinates, respectively. 

^We are ignoring the fact that globally U{Nv + N\Nv) = [SU{Nv + N\Nv) (g) U{1)]/Zn, i.e. it is a 
coset rather than a direct product. This is irrelevant for small transformations. 
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As discussed below, the difference between this group and eq. (|T0]) is the mathematical 
result which underlies the need to include the $o field in the quenched, but not the PQ, 
chiral Lagrangian. 



3 Symmetry breaking and the need for <l>o 

As could already be seen in the previous section, the introduction of valence and ghost 
quarks modifies the standard flavor symmetry structure of the theory, and consequently 
the formulation of the chiral effective theory. One striking example is the non-decoupling 
of $0 in the fully quenched theory. The main goal of this section is to analyze this 
phenomenon in some detail, and to demonstrate why it does not carry over to PQ theories. 
To study this, we investigate which two-point correlation functions contain poles at low 
energies. In other words, we find the degrees of freedom that one must include in a low 
energy effective theory for quenched and PQ QCD. 

Throughout this section we consider the PQ theory in the chiral limit. This allows 
us to use the condensate as an order parameter, and to use Goldstone's theorem to 
determine which channels have massless poles. In the usual way, the chiral limit is to 
be approached by working at non-zero quark masses and then taking the masses to zero 
after the volume has been sent to infinity. In the PQ theory there is, however, a subtlety 
concerning the chiral limit. As noted in Ref. |^, chiral perturbation theory predicts that 
the PQ theory is singular if one sends valence quark masses to zero with fixed non-zero 
sea-quark masses. In particular, the condensate itself is singular in this limit. This 
singularity is similar to that which occurs in the quenched theory |T3|, [15|. To avoid 
this singularity, one should send all the quark masses to zero simultaneously with fixed 
ratios.0 

Such a limiting proceedure is not available for the quenched theory. Thus, in the 
following, when we refer to the quenched theory, we have in mind working close to, but 
not in, the chiral limit. Since our focus here is on the PQ theory, we do not revisit the 
subtleties associated with the quenched chiral limit. 



3.1 Symmetry breaking pattern 

As in QCD, we choose the vacuum expectation value (VEV) 

^ab = (QaQb) (12) 

as an order parameter for chiral spontaneous symmetry breaking. 

We first consider the vector symmetries. It was shown long ago by Vafa and Witten 



that vector symmetries do not spontaneously break in vector-like gauge theories ||16 
Their derivation does not make use of Hilbert space states and operators, and relies 
only on the fact that the quark determinant leads to a real and positive measure in 
the functional integral over the gauge fields. This is still true for quenched and PQ 



This limiting proceedure also avoids divergences from exact zero modes in topologically non-trivial 
sectors. 
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QCD, and the Vafa-Witten result still holds. Consequently, ^lab is invariant under vector 
transformations. It is in fact easier to see what this implies for the related quantity 

nab = (QbrQar) (13) 

(r is a Dirac-color index that must still be contracted to form a Euclidean scalar). This 
transforms under ( "fake" ) vector transformations in the following way: 

n^Vnv\ V e SU {Nv + N\Nv) . (14) 



The invariance of under eq. (Ill) leads to 

Q = Uj6ab (15) 

where u; is a constant. Interchanging Q and Q fields, we obtain 

Qab = -UJ^abSa (16) 

where we introduce the notation 



1 a is a quark index, 
-1 a is a ghost index. 



:i7) 



The question of whether the axial symmetry is spontaneously broken or not depends 
on the value of u. The sea sector of PQ theories is equivalent to unquenched QCD with 
quark flavors. An important implication of this fact is that the spontaneous breakdown 
of axial symmetries in QCD, signaled by the non- vanishing of (qq), is duplicated in the 
sea sector of PQ theories, thus implying that u ^ 0. In other words, given that there is 
spontaneous chiral symmetry breaking in QCD, the symmetry breaking pattern of PQ 
QCD is known:g 

SU{Nv + N\Nv)l ® SU{Nv + N\Nv)r ® f/(l)y ^ SU{Nv + N\Nv)v ® U{l)v (18) 

This argument does not carry over to the quenched theories, because there is no QCD- 
like sea sector. For quenched theories the spontaneous breaking of axial symmetries is 
an additional assumption — though one that is supported by numerical evidence. If we 
make that assumption, the symmetry breaking pattern for quenched QCD is 

[SU{Nv\Nv)l ® SU{Nv\Nv)r] k U{1)v ^ SU{Nv\Nv)v x U{l)v. (19) 



3.2 Low energy poles from symmetries 

Once the symmetry breaking pattern has been established, it is a standard result in field 
theory that the number of pseudo-Goldstone bosons is determined by the number of non- 
anomalous generators that act non-trivially on the vacuum. For the case at hand one thus 
expects (2iVy + A^)^ — 1 Goldstone particles for PQ theory. Though the group structure in 

^Here we are assuming N > 2, since there is no chiral symmetry to break for = 1 QCD. 
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the quenched case is shghtly different, the counting argument still gives {2Nv)'^ — 1, which 
is the same expression as for PQ theories with = 0. There is, however, a significant 
difference between the two cases. As the following more careful analysis shows, the simple 
counting argument correctly predicts the number of Goldstone particles only in the PQ 
case. For the quenched case it turns out that there are (2iVy)^ fields which exhibit long 
range correlations. The additional field that is needed for quenched QCD is none other 
than $0- 

Consider the two-point correlation function between an axial current (x) (defined 
in eq. (|115|) ) and a pseudoscalar density 0'^"^'^(O) = Q75T'Q(0), 



^r'Hx) = (jf;(a:)0(^')(O)). (20) 

Here T and T' label generators of the symmetry group. From Euclidean invariance, its 
Fourier transform must have the form 

C^^''''\p) = ^P,F^^'^'\p'), (21) 

where F is an unknown function. Next, consider the Ward identity which follows from 
applying an infinitesimal axial transformation with generator T to (0'-"^ ^(0)): 

d,Cl^'^'\x) = -5{x){5T<P^^'\0)). (22) 

(for the definition of S]^ (p see eq. ( |116|) .) Together, eqs. (^T|-^2|) give 

p^Fm(p^) = -(5jV(^')(0)). (23) 

A straightforward calculation of the right hand side, with the use of the VEV (0), gives 

/F(^'^')(/) = 2a;str(TT'). (24) 

(Note that the order of T and T' is important for graded generators.) Eq. (2^) relies on the 
fact that T generates a true (non-anomalous) symmetry of the theory. For quenched and 
PQ theories the non-anomalous symmetries satisfy str(T) = . Consequently, ' (p) 
has a pole at p = for any T and T' for which str(T) = and str(TT') ^ 0. 

It is straightforward to show that, for both quenched and PQ theories, all the flavor 
off-diagonal ("charged") generators give rise to light poles, like in QCD, and that they do 
not mix with the neutral ones. Subtleties only arise in the neutral sector, and to study 
them we use a specific choice of diagonal generators that is presented in [Appendix B| . 
There we find that, for the PQ theory. 



str(T,rb) = diag(l, . . . , 1, -1, . . . , -1, -1, 1). 




Here the Ta are diagonal members of the algebra, all straceless but the last one, T2NV+N 
/, which generates the anomalous U(1)a- Applying eq. ( |130| ) to eq. ( p^ shows that a 
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complete set of (flavor neutral) fields that give rise to long range correlations in two-point 
functions is 



2Nv + N - 1. 



(25) 



Note that for A'^v" of these fields, the coefficient of the massless pole has an unphysical 
sign, corresponding to the entries with —1 in eq. ( p.3CI| ). The symmetries do not imply 
that (j)^'^2Nv+N) — has a pole at p = 0. Consequently, to describe the long 

distance parts of the two-point functions in the PQ theory we do not need to consider 
correlators involving the corresponding field, $o- 
In the quenched case eq. ( |130| ) is replaced by: 



/I 



str (TaTft) 



V 



1 

1 0/ 



Again, the generator of the anomalous U(1)a is the last one, T2Ny oc / (defined in 
Appendix B ). All the other generators are straceless, including T2NY-1 = I- obtains 
a new non-trivial equation from eq. (p^ by choosing T = I, which is non-anomalous (a 
fact unique to quenched QCD), and T' = I. We thus find that the Fourier transform of 



has a pole at p = 0. Since $0 has the same quantum numbers as 



(26) 



(27) 



we conclude that it must be included in the effective Lagrangian in order to correctly 
reproduce the low energy behavior of quenched QCD. 

Let us summarize why one needs to include $0 in quenched, but not PQ, chiral per- 
turbation theory. A field should be included if either the corresponding pseudoscalar 
density or axial current can be shown to couple to a Goldstone boson. In unquenched 
and PQ theories, the argument given above shows that the coupling of a density with 
flavor generator T implies the coupling of the corresponding current [since sti {TaTf,) is 
diagonal]. The peculiarity of the quenched theory is that the current with the same quan- 
tum numbers as $0, i-e. the flavor singlet j^^^, is not anomalous and thus should couple 
to a Goldstone boson. This is despite the fact that the corresponding density, Q^i^Q, 
is invariant under non- anomalous transformations. Conversely, the anomalous current, 
j^^^, is not a flavor singlet, and so, as shown above, the corresponding density couples 
to a Goldstone boson, contrary to naive expectation. Mathematically, this peculiarity 
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is due to the fact that str(J) = 0, which leads both to the key result that str (TaT;,) is 
not diagonal, and to the semi-direct product structure of the quenched symmetry group, 



eq. (11). This is the sense in which the symmetry structure of the quenched theory leads 
to the need to include the $o- 

In [Appendix C| we give an alternative argument for the need to keep $o in quenched, 
but not in PQ, theories. 



4 Effective Lagrangian 

We now turn to the construction of the low energy effective Lagrangian for PQ QCD. We 
first recapitulate the standard approach, based on the "fake" symmetry group discussed 
previously. We point out the problems with this approach, and spend most of this section 
discussing the extent to which they can be alleviated. Our conclusion is that the standard 
approach is appropriate when doing perturbative calculations, although the justification 
for using the effective theory is considerably less rigorous than for QCD itself. Some 
of these points have been discussed previously in the random matrix literature, e.g. in 



Ref. [|Tl|], although in the context of effective theories including $o- Our aims here are 
to clarify this work, and to generalize it to the case at hand in which $o is absent. Parts 
of our discussion are based on an analysis of the quenched effective Lagrangian given in 



Ref. 17. 



4.1 The standard approach 

The fields in the effective theory can be limited to those describing the pseudo-Goldstone 
hadrons. These, as we have seen, are in one-to-one correspondence with the generators 
of 

[SU{Nv + N\Nv)l ® SU{Nv + N\Nv)r]/ SU{Nv + N\Nv)v • (28) 



In the standard approach [|I0| the pseudo-Goldstone fields lie in this coset space, with 
the parameterization 

^w--(^)' IS)' <^^' 

Here the blocks correspond to quarks and ghosts, and rji and 772 are independent Grass- 
mann matrix fields. We are free to interpret 772 as in which case U is unitary. The 
absence of $0 implies the constraint str($) = ticj) — ticj) = 0. The symmetries are imple- 
mented as usual: 

U LUR^ , (L, R) G SU{Nv + N\Nv) , (30) 

with the VEV (U) = I breaking the symmetry in the required way (note that U trans- 
forms like Q — see eqs. (|13|-|1^)). The invariant effective Lagrangian is 

C {U) = y str [dUdU^) - y str [xU^ + Ux) + higher orders, (31) 
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where x is proportional to the quark mass matrix, x = 2BoM, and Bo and / are unknown 
parameters. Clearly this development mirrors that for QCD step by step. (Note that 
eq. ( pID differs from the form of the chiral Lagrangian for PQ QCD which is usually 
quoted: The absence of $o nieans that arbitrary functions of $o which might multiply 
each term in the PQ chiral Lagrangian are absent.) 

A problem with this effective Lagrangian becomes apparent when developing pertur- 
bation theory by expanding U about its VEV. The "str" implies that the fields have 
kinetic and mass terms with the wrong sign, so that (U) = J is not a minimum of the 
action.0 This is dealt with in the standard treatment by simply ignoring the instability. 
A justification for the use of "wrong sign" propagators for some of the mesons is that it 
implements cancelations which correspond at the quark level to the desired cancelations 
between valence quark and ghost loops. Actually, once $o is excluded, the connection be- 
tween pseudo-Goldstone propagators and underlying quark and ghost flows is less clear, 
so this qualitative justification is less convincing. Clearly a better-founded treatment is 
desirable. 

A related concern with the standard approach is that the chiral symmetry group of 
the PQ QCD functional integral is not SU{Nv + N\Nv)l ® SU{Nv + N\Nv)r, and the 
unbroken vector subgroup is not SU{Nv + N\Nv)- As explained in [Appendix A| , the 



full group is the subgroup of SL{Nv + N\Nv)l ® SL{Nv + N\Nv)r, in which (using the 
nomenclature of Ref. |]l8l) the "bodies" (non-nilpotent parts) of the left- and right-handed 



transformations are related by 



L = exp(^$i), i? = exp(^<l>^), {'^ l) gg\l,^^ = {'^ r) ggi,^^ , (32) 



(with str ($L,/?) = 0) while (as explained at the end of [Appendix A| ) the vector group is 



the subgroup of SL{Nv + A^|A^y) in which the body of the ghost-ghost part is unitary: 

L = i? = V = exp(^$y) , i^v)gg\l^^ = {^v)ggL^^. (33) 

The coset of these symmetries can be parameterized by transformations with L = R~^, 
so that 

L = R-^ = A = expt^A , (*A)33|Ly = -('^'A),,|body • (34) 

Note that aside from the condition on the body of the ghost-ghost block, is an 
arbitrary, straceless matrix. One would then expect that the correct Goldstone fields live 
in this coset space, i.e. 



, str(<l>') =tr(0') -^tr(0') = 0. 

body 



(35) 



^Strictly speaking, this comment applies to fields living entirely within (j). However, the anomaly 
constraint implies that one neutral meson must have components from both (f> and (f) — that corresponding 



to the generator T2NY+N-1 in Appendix B. Since str(r2Ar -+n-i) ^ ^^^^ mixed field also turns out 



to have a kinetic term with the wrong sign. 
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(Note the crucial factor of i multiplying (p'.) cj)' is an arbitrary matrix of complex c- number 
fields, constrained only by the stracelessness condition on 

It turns out that neither eq. ( p9D nor eq. (|35D is correct, but before discussing this 
point it is useful to understand how the form of the effective Lagrangian would differ were 
the Goldstone fields given by eq. (^). This field transforms like U' LU'R~^, with 
L, R given by eq. (0). The invariant Lagrangian is constructed from U', U'~^ , the mass 
term (a spurion which transforms like U'~^), and other sources. The rules for combining 
these are in one-to-one correspondence with those for the usual chiral Lagrangian (with 
U ^ U\ W ^ U'-^). Thus one finds that the most general Lagrangian has the standard 
form, (pID, except that U'~'^ appears rather than . We stress that this holds to all 
orders in the chiral expansion. The fact that sdetf/' = 1 implies that there are no 
additional $o-like terms. One might be concerned that the resulting Lagrangian has, in 
general, an imaginary part. This will turn out not to be true in the final form we consider 
and so we do not discuss it further here. 



4.2 The correct Goldstone manifold 

Now we return to the issue of the correct Goldstone manifold to use. Clearly it should 
be based on the true symmetries, and thus contained within eq. (|35|). The problem is 
that U' in eq. (|35D has too many degrees of freedom, and we must pick an appropriate 
sub-manifold. This issue has been addressed by Verbaarschot and collaborators, in the 
context of a theory which contains $o (see e.g. Ref. [0). They argue, based on the 
mathematical results of Zirnbauer ||T^, that the appropriate integration domain for the 
effective theory is the maximal super-Riemannian manifold contained in the coset of the 
true symmetry group and the unbroken vector group. This results in a Goldstone field 
parameterized hyf] 



U"{x) = exp ( j-^j , 



\r]'^{x) i4>"{x)J ' (36) 



This form is then inserted in the effective Lagrangian described in the previous paragraph, 
i.e. (|3lD with U U" and — >■ U"~^ . This leads to a Lagrangian whose body is real. 



In the following, we discuss the reasoning leading to the form (|36|), and address a 
technical difficulty which arises when extending the argument to the theory without $o- 
We then discuss the theoretical foundations of PQ chiral perturbation theory, and along 
the way give further physical motivation for the choice (pG]). 

The important differences between the parameterizations (^), (P3| ) and ( PUj) are in 
the ghost-ghost and quark-quark blocks of the exponents.^ We begin by focusing on the 

^For explicit non-perturbative calculations, the authors of Ref. | [Tl[ | use a different parameterization. 
We use the form (|3^) as it is closer to the standard parameterization of eq. ([29|). 

®The anticommuting parts of the exponents have the same form in all three parameterizations. We 
note that the constraint 771 = ri\ imposed in ( p9| ) does not change the rules for Grassmann integration, 
in which rji and 772 are treated as matrices composed of independent Grassmann variables. There are 
also no issues of convergence in the Grassmann integrals. 
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former. Here the difference between tlie full coset space of eq. (|35| ) and the manifold in 
eq. (|36| ) is only that the soul of (f)' is not required to be hermitian while that of 0" is. Thus, 
0' appears to contain more degrees of freedom than 0". However, since Gaussian integrals 
over c-numbers are independent of the soul, as long as the integrals are convergent (see 
Appendix A| ), these extra parameters can be absorbed into 0'. Thus the two forms ( |35[ ) 



and ( p6D do not differ in this respect. 

They do, however, both differ from the standard parameterization, eq. (p9D , by the 
extra factor of i. This factor resolves the stability problem raised above. If we expand 
either U' or U" around / to quadratic order, then the extra factor of cancels the 
minus sign from the supertrace, and the kinetic terms for 0' or 0" (which are equivalent 
parameterizations at this order) have the correct sign for stability. This also fixes the 
overall sign of the kinetic term in the Lagrangian. The same discussion holds for the 
mass term, as long as M is proportional to the identity. We postpone discussion of the 
(most relevant) case of non-degenerate masses until later, since to address this we need 
first to understand the effect of the absence of $o- 

We now turn to the quark-quark block of the exponents. The full coset space 



contains extra fields in this block, compared to the standard form of (p9D , since 0' is 



not constrained to be hermitian. These extra fields correspond to the generators of the 



non-compact part of the broken symmetry group. The choice of Ref. [|TT| is to keep only 
the usual, hermitian part of 0'. A technical reason for doing so is that, if one expands 
about U' = I, the kinetic term in the action (the sign of which is now fixed) is minimized 
in these "hermitian" directions, but not in the "non-compact" directions. Thus, with 
this choice, the propagators of both 0" and 0" have the correct signs. For the moment 
we will accept this as sufficient reason for making this choice, but return to the point 
below when we discuss the foundations of PQ chiral perturbation theory. 

4.3 Constraints from absence of $o 

The analysis so far has not accounted for the constraints imposed by the absence of $0; 
by the fact that only non-anomalous symmetries should be represented by the effective 
Lagrangian. There is no problem when we consider the full coset space of eq. (|35|) . The 
constraint is that str($') = 0, and this can be satisfied by expanding the diagonal or 
"neutral" part of $' using the basis for diagonal straceless generators given in Appendix 

2Nv+N-l 

Keu= E (37) 
a=l 

The important point for us is that while all but one of these generators are contained 
entirely within either the quark or the ghost sectors, there is one generator (T2Nv+n-i 
in our ordering) with components in both sectors. Now consider the restricted manifold 
of eq. (|36D. Here the constraint str($") = tr0" — 2tr0" = is made more stringent by 
the fact that 0" and 0" are hermitian. Thus 0" and 0" must be separately traceless, and 
there can be no component of $'^g„ proportional to T2Nv+n-i- Thus the restrictions of 
the manifold remove not only $0) but also another neutral generator. This is a problem 
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because, as seen in the previous section, there are long range correlations in channels 
with the quantum numbers of this generator. In other words, U" is missing one of the 
neutral pseudo-Goldstone bosons. 

This problem can be solved by a small change in U" . As we have seen, U" was 
constructed so that all bosonic fields have correct sign propagators. For the neutral 
fields, this requires taking components in the quark sector to be real, but those in the 
ghost sector to be imaginary. The generator T2NV+N-1 is problematic because it straddles 
both sectors. However, since it satisfies str(T^) = — 1, and is thus ghost-like, we propose 
that the corresponding field, ^"^A^y+Af-i^ should also be taken imaginary. Explicitly, we 



propose replacing in eq. (|36D with 

c/i ' ne 



Nv+N-1 2Nv+N~l 
0=1 a=Nv+N 



(38) 



Here "c/i" refers to the off-diagonal, or charged, parts of the field. In this way we include 
all the neutral Goldstone fiavors, while maintaining the anomaly constraint, and the 
condition that all kinetic terms have the correct sign. 

This is not quite the end of the story. The convergence of the functional integral 
depends also on the mass term. Given the choice of Goldstone fields just described, the 
Goldstone boson mass matrix can be calculated, and the functional integral converges 
only if its real part is positive definite. As one might expect, there is competition between 
the sea quark masses and the valence (and ghost) ones. It turns out (the derivation is 
somewhat tedious but straightforward, and we do not include the details here) that the 
quark masses must satisfy 



NvXviNvW + Nx^) < {Nv + Nf. (39) 



Here xv and xs are the average valence and sea quark masses, and Xv^ the average 
inverse valence quark mass. Restricting lattice simulations to satisfy this constraint is 
undesirable and very likely unnecessary. Although the theoretical description of the 
simulations in terms of the PQ chiral effective theory is ill defined when relation eq. ( |5^ 
is violated, there is nothing special about this point in parameter space as far as the 
underlying PQ QCD is concerned. Thus, the simulations will show no special behavior 
when eq. (^) is violated. In addition, chiral perturbation theory is insensitive to these 
global convergence considerations, and if it is regular at this point, and provides an 
adequate description of low energy PQ QCD when the quark masses are such that the 
inequality is satisfied, it is unlikely that the chiral expressions derived from this theory 
will suddenly cease to make sense just when the inequality becomes an equality. 

Aside from this point, our construction has yielded an effective Lagrangian with the 
correct number of Goldstone fields (i.e. the same number as there are independent Ward 
identities), and which one can consistently expand about {U") = I. While in principle 
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one can work with this Lagrangian, it is more convenient for perturbative calculations to 
change variables as follows: 



// 

ch 



a 



-ia", a = Nv + N,2Nv + N -1, 



(40) 



a", a = l,Nv + N-l, 



The effect of this change is that U" now appears to have the form of U in eq. (^91) . Since 
the underlying integrations are unaltered, the only effect of this change of variables is to 
shuffle factors of i between propagators and vertices (and possibly external fields). For 
example, all the charged (p propagators now have an overall "wrong" sign, while vertex 
factors appear to satisfy the graded SU{Nv + N\Nv) symmetry. Indeed, we can now 
pretend that U G SU^Ny + N\Ny), and ignore the subtleties of this section, as long as 
we use the prescription that wrong sign kinetic terms lead to wrong sign propagators. 

4.4 Foundations of PQ chiral perturbation theory 

As promised, we now return to the foundations of PQ chiral perturbation theory. We 
first revisit the issue of the extra "non-compact" generators in the quark sector. Since 
this sector is a simple generalization of QCD, the construction and justification of its 
effective low energy theory is well understood (see, e.g. Refs. [||, ^). In particular, the 
issue of extra generators arises also in QCD itself, though it is rarely discussed. 

It is useful to recall the steps needed to construct the effective Lagrangian for QCD- 
like theories with Nf flavors. First, one derives the Ward identities following from the 
chiral symmetries. Second, given the assumed VEV, Goldstone's theorem determines that 
there is one light pion field for each axial current. Third, one assumes pion dominance of 
correlation functions, i.e. that the light fields are the only relevant degrees of freedom. 
One then writes down the most general local Lagrangian incorporating these degrees of 
freedom which is invariant under the chiral symmetries, and asserts that this will yield 
the most general amplitudes consistent with the Ward identities and the principles of 
relativistic quantum mechanics . To formalize this, one introduces sources for currents, 
scalar and pseudoscalar densities, etc., and defines the generating functional in the usual 
way 



The Ward identities of the theory are encapsulated in the invariance of Z under gauge-like 
transformations. 



A similar definition is used for Zes[f^,f^,S,P,...]. The claim is that, by adjusting 
the coefficients in the effective Lagrangian, including contact terms, one can match the 
generating functional between QCD and the effective theory. 




(41) 



/i: ^ L/^L-i - t[d^L)]L-' , f>^ ^ Rf^R-' - t[d'^R)]R-' , 
S -iP^ L{x)[S - iP]R-\x) , S + iP R{x)[S + iP]L-\x) . 



(42) 



Zeslfl fi^,S,P,...]= Zqcd[/^, fj^,S,P,...], 



(43) 
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to any desired accuracy in a momentum and quark mass expansion pO |. 

The issue at hand is what group the transformations L and R belong to. This de- 
pends on whether Zqcd is defined using the operator formulation or with a functional 
integral. In the former case, the transformations of ip and tp'^ are related, since ip = ip'^'jo. 
Thus only unitary chiral transformations are symmetries [as in eq. (^], and L{x) and 
R{x) in (^2]) are elements of SU{Nf). On the other hand, in the functional integral 
formulation, the fact that ip and ip are independent variables leads to a larger symmetry 
(see [Appendix A|) : L,R E SL{Nf). Which of these symmetries should Z^s respect? The 
answer is the smaller, unitary symmetry. This can be seen in two ways. First, the "extra" 
transformations contained in SL{Nf)/ SU{Nf) lead to non-hermitian sources S and P 
[see eq. (^2])], and so to a non-hermitian Hamiltonian. Thus they move us out of the space 
of physical theories. Second, as shown in [Appendix A| , the extra transformations do not 
lead to additional Ward identities, and thus do not lead to additional Goldstone bosons. 
Since the effective Lagrangian should certainly respect the symmetries of the operator 
formulation of the theory, and this leads to all the desired Goldstone bosons, one should 
not extend the effective theory to include non-compact symmetries. We conclude that 
the "normal" choice of pion fields in the quark block of the Goldstone matrix, eq. (p6[), 
is correct. The only exception is for the part of 0" proportional to the identity in this 
block (i.e. the field (j"2Afv+JV-i^_ rpj^^ argument just given does not apply to this field 
since it is a fiavor singlet in the quark block. 

Finally, we address the extent to which one can derive the effective Lagrangian in 
the graded sector of the theory (i.e. the parts involving ghosts). The first two steps 
followed in the quark sector go through in this sector as well: the Ward identities of 
PQ QCD are the graded generalizations of those in QCD, and the symmetry breaking 
pattern in PQ QCD follows from that in QCD. From these results, we have established 
the presence of massless poles in two-point functions of broken symmetry generators. 
In the standard approach, the next step is to interpret each pole as being due to a 
physical single particle state created by the corresponding operator, from which follows 
the result that any correlation function including these operators will have poles at the 
same positions. This deduction is key in justifying the effective chiral Lagrangian for 
QCD, but its extension to PQ QCD is not obvious. Furthermore, when one includes 
quark masses as small perturbations in QCD, one concludes from similar arguments that 
this results in a small shift in the position of the poles. In PQ QCD, on the other hand, 
quark masses have a substantially different effect-they can lead to the appearance of a 
double pole at = with a coefficient proportional to the quark masses, as will be seen 
explicitly in sec. ^. 

What is lacking in the PQ case is a positive-norm Hilbert space interpretation, since 
the theory contains ghosts. This is not to say that a derivation of the effective Lagrangian 
is not possible for PQ theories, but rather that a generalization of the standard methods 
is needed. Since quark and ghost correlation functions differ simply by signs, we speculate 
that the appropriate theoretical framework only differs from the standard one by requir- 
ing a vector space with indefinite metric (e.g. single anti-ghost states having negative 
norm), as well as bosonic ghost creation and annihilation operators. The Hamiltonian 
would still be physical — having real eigenvalues bounded from below — and ghost opera- 
tor commutators would be causal. In this way Lorentz invariance is maintained, while 
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unitarity is lost because of transitions between positive and negative norm states.^ It is 
perhaps plausible in such a set-up that a generalization of standard "pologogy" could be 
derived, in turn leading to an effective Lagrangian, and that the requirement of Lorentz 
invariance would force this Lagrangian to be local. 



5 Integrating out the <l>o 

The work of the previous section has shown that we can write the effective chiral La- 
grangian for the PQ theory in terms of an SU{Nv + N\Nv) field U, transforming as in 
eq. (|30|). In this section we need only write the form of this Lagrangian schematically: 



C{U) = J2WU), (44) 

i 

where the ii are unknown parameters, and Oj are operators constructed from U, W, 
their derivatives, the rescaled mass matrix and sources for external operators. We 
note that the operators allowed by the graded chiral symmetry are identical in form to 
those allowed by the usual chiral symmetry in QCD. In particular, there are no additional 
operators. 

An important property of PQ chiral perturbation theory is that correlation functions 
involving external sources restricted to lie within the sea-quark sector are identical to 
those obtained using the effective chiral Lagrangian for unquenched QCD (with Nf = N, 
and without the rj'). This identity, which is trivial at the quark level, can be seen 
diagram by diagram in chiral perturbation theory, due to cancelations between diagrams 
in which a valence quark is replaced with a ghost quark. A compelling general argument 
in support of this fact, though not a proof, is the following. Let Ai be the mass matrix 
for a specific valence quark and the corresponding ghost, A4 = diag(my, my). Let A be 
an operator which does not depend on this valence-ghost pair. We write its expectation 
value as (A) = a{Ai), making only the Ai dependence explicit. The flavor symmetry 
group includes SU(1|1) transformations that mix the valence and ghost quark fields. 
These transformations are equivalent to a change of variables in the functional integral, 
in addition to a group transformation (by conjugation) of Ai. Since A is unchanged 
by the change of variables, it follows that a{A4) can depend only on SU(1|1) invariants 
constructed from Ai. These, however, can only be supertraces of powers of the mass 
matrix str(A^") = 0, so a is independent of the valence (and ghost) quark mass. To show 
independence of a on the existence of the valence-ghost fields in the theory one still needs 
to show that any effect from these fields introduces A4 dependence. 

It follows that, with external sea-sector sources, one can restrict the internal U field 
to take the following block-diagonal form 

Uqcd = {In^,SU{N),Im,) . (45) 



^This is indeed the structure observed in explicit calculations of scattering amplitudes in quenched 
chiral perturbation theory: the results are Lorentz invariant but not unitary pl|. 
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Since the terms in C{U) are in one-to-one correspondence with those in the QCD effective 
Lagrangian, and the restriction (^51) does not change the form of the operatorsj^ it follows 
that both PQ and unquenched Lagrangians share the parameters As we emphasized 
in H, this means that the parameters ii, which describe the chiral expansion for 
the unphysical PQ theory, are in fact physical. For example, at NLO, ii are the GL 
coefficients which encode our current experimental knowledge (and ignorance) of QCD 
at low energies for the light mesons. We stress that the ii do depend on N, i.e. the GL 
parameters that one obtains depend on the number of light sea quarks. 

From a theoretical point of view, this completes the construction of the PQ chiral 
effective Lagrangian, and the demonstration that it contains only physical parameters. 
When doing perturbative calculations with this Lagrangian one must, however, imple- 
ment the constraint that sdet U = 1, or, if we write U = exp(2z$//), that str($) = 0. 
A standard way of achieving this is using straceless generators. For the realistic case 
of Nf = 3 and considering only mesons, A'^v' = 2, calculations involve the generators of 
SU{5\2), and can be quite tedious. In this section we point out that a simple alternative 
is to reintroduce the $o field — as a calculational device and not as a physical field — and 
then to integrate it out. This allows one to work in the "quark basis", rather than with 
the actual pseudo-Goldstone fields. Calculations in the quark basis are more transparent 
since one can trace quark fiow through each diagram, and see the cancelations between 
valence quarks and ghosts very simply 0, Keeping $o also allows us to reinterpret 



previous calculations, which have included it as a physical field, as applying to the theory 
without this field (a point discussed at the end of this section). 

5.1 Functional integral approach 

We reintroduce $0 by enlarging [/ to S G U{Nv + N\Nv), 

E{x) = U{x)exp(^^^) , (47) 



fVN 

and considering a theory with Lagrangian 

£'(E) = ml^l + CiJ:) (48) 
= ml^l + J2^^0^i^)■ (49) 



-"With eq. ([45|), %, U and dU are proportional to matrices of the general form 

(46) 




where the three entries correspond to valence, sea, and ghosts, and A, B stand for any Ny x Ny and 
N X N blocks, respectively. Any product of such matrices still has this structure, which causes each 
supertrace in the Lagrangian to reduce to a simple trace over the sea-sea block. 
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In words, we simply replace U with S in all the terms in C{U) and add a mass term for 
$0- It is useful to make the $o dependence explicit by expanding S: 

oo 

£(S) = C{U) + ^ Rj{U, d){%y , (50) 
i=i 

where we allow Rj{U,d) to contain derivatives acting on the $o field. The enlarged 
theory is then 



DSexp l^-y /:'(S)J (51) 
DUD^expl^- J (^CiU) +ml^l + Y^RAU){^oy'^y (52) 
We assume that the theory is regulated with a chirally invariant fixed cut-off, such as 



the lattice [23 



We now take the limit mg oo, and argue that we then return to the original theory 
with E — i> [/, i.e. the theory that we want to do calculations with. The argument goes as 
follows. We expect the fluctuations in $o to be of (9(l/mo), and thus that all the Ri terms, 
which are independent of mo, are suppressed compared to the mass term. Similarly, 
correlation functions calculated in this theory should become independent of Ri at mo 
oo. Furthermore, the expectation value of any operator 0(S) = 0{U) + X]j '^i(f^)(^o)"' 
satisfles 



{Oi^))c' = j ^^^^^ ^(^) (- / ^'(^)) (53) 
^ -y^jDUD%0{U)e^^[-!{C{U)+ml^l)\ (54) 

= ^jDUO{U)ew{-jc{U)^ (55) 

= {0{U))c. (56) 

So we obtain correlation functions in the theory we want, with any $o contribution to 
the external operators being projected out. 

It is crucial for this argument that the theory be regulated in such a way that loop 
momenta are limited.^ This allows $o to be integrated out in a trivial way. Without a 
flxed cut-off, the (9$o)^ term [implicitly contained in i?2(f^)^o] dominate over the 
mo term in loop integrals. This in turn leads to a non- decoupling of $o- For example, 
using dimensional regularization, $o tadpole diagrams give contributions proportional to 
(mo//^) ln(mo//iDR). With a flxed cut-off, by contrast, the contribution is of the form 
A^/(/^mo), and vanishes when mo — > oo. 

There is a second subtlety which could invalidate the argument just presented. In 
physical theories heavy particles "decouple" , meaning that physics at energies much lower 
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than their mass is independent of the details of their interactions and dynamics. This 
idea is at the core of effective field theories in which only light particles are included as 
explicit dynamical degrees of freedom. In this sense, taking the mass of a particle to 
infinity in physical theories is well defined. For unphysical theories, however, one might 
be concerned that this limit does not exist. That this is a legitimate concern is shown 
by the fact, explained below, that the limit cannot be taken in quenched QCD. We will 
see, however, that this is not a problem for the PQ theory. 



5.2 $0 in the flavor- neutral propagator 

To understand this point, and also to gain insight into the nature of the ttiq — ^ oo limit, 
we consider the form of the propagator of "neutral" mesons, i.e. those created by the 
diagonal elements of $. It is sufficient to consider this propagator since this is the only 
place that the ttIq term enters when one develops perturbation theory. We begin with 
some notation. We use Il{x) to denote the pion field including <I>o: 

E = exp (^^^ , n(a;) = + M^)^ ■ (57) 
The set of meson fields {Tfij} that make up the quark basis is defined through 

2Nv+N 

n(x)= j2 T^J7l''i^)^ m,)ki = 5,k5ji, (58) 
and the flavor neutral part is therefore 

2Nv+N 

n„eu(^)= Yl Tu7l'\x). (59) 
i=l 

The neutral part of 11 can also be decomposed using the basis of generators of U {Ny + 



N\Nv) given in [Appendix B 



2Nv+N 

nneu(x)= J2 Ta(y\x), (60) 
a=l 

where the generators satisfy 

str (T„Tb) = Qah- (61) 

The first 2Nv + — 1 components of a" are true pseudo-Goldstone fields, while the last 
entry is the field added by hand, $o- 

We consider all the quadratic terms in the action that contribute to the leading order 
propagator. At this point we do not specify what these terms are, only require that the 
$0 mass term is included. For a given momentum, the inverse propagator will have the 
following form: 

G7\ = 2^ }\. (62) 
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Here we have separated off the last row and column of the matrix, corresponding to 
entries involving $o- Thus A [a {2Ny + N — 1) x {2Ny + — 1) matrix] is the inverse 
propagator in the theory without $o- The important point is that A, B, C and d are 
independent of mo, and finite (because of the momentum cut-off). Using this, one can 
easily show ( [Appendix D| ) that the only requirement needed for the propagator to have 



a limit as itlq — oo is that A be non-singular. The limit is then 



— 'A'^ 0^ 

Gw) = lim G'(^) = ( ^ ^ ) . (63) 



This propagator has the two properties needed so that one is effectively doing the calcu- 
lation using only the physical SU{Nv + N\Nv) degrees of freedom: (1) its projective form 
removes $o from the theory — factors of $o in vertices and external fields simply do not 
propagate. (2) the propagator in the physical subspace, A~^, is the correct propagator 
for these degrees of freedom alone. |^ 

Thus the only remaining question is whether A~^ exists. To study this we first 
consider only the kinetic + ttIq terms in the Lagrangian eq. (^9]) . A can be easily read 
off, using 

str (dUneudUneu) oc Qabda^da^ , (65) 
and the explicit form of gab from [Appendix B|: 



PQ : A = / diag(l, . . . , 1, -1, . . . , -1, -1) 
Q: A = /diag(l,... ,1,-1,... ,-1,0) ^ 

Clearly A is invertible only in the PQ case, and the limit mo oo can be taken only in 
that theory, but not in quenched QCD. 

The usual mass term in the Lagrangian, which is of the same order in chiral pertur- 
bation theory as the kinetic term, can be treated as a vertex. This leads to a geometric 
series of tree diagrams, all contributing to the leading order. The propagator in these 
diagrams, when it exists in the mo oo limit, removes any $0 contributions. This 
means that the conclusions derived above from eqs. ([66|) still hold. There is, however, 
a potential loophole in this argument: the infinite sum over tree level diagrams and the 



One might be concerned that the fact that some of the fields have the wrong metric might lead to 
a basis dependence in the results of calculations. Indeed, the transformation between the "sigma" and 
"pi" bases, a" — X^f, ^''ftTr'''', is not unitary, but instead is a generalized Lorentz rotation. This follows 
from the results 

str (nL J = """dabo' ^ E ^''''~9ab^'' ' (64) 



where gab (defined in Appendix B) plays the role of the metric tensor, and gab has the same signature 
but with elements permuted. Since Lorentz-like rotations do not preserve Euclidean inner products, 
projections onto subspaccs arc not in general invariant under these rotations. This turns out, however, 
not to be relevant in our calculations. All that we need in order to show the basis independence of our 
calculations is that the propagator of the matrix field, (n„eM(a;)n„e„(0)), which is the building block of 
perturbation theory, is the same in either basis. It is straightforward to show that for this to hold, it is 
sufficient for A to be invertible, which is clearly the case. 
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mo oo limit might not commute. That this is not a problem can be seen directly from 
previous calculations of the LO propagator in chiral perturbation theory in which the 
quark mass term was included before taking the mo — > oo limit. In the case of PQ QCD 



with 3 sea quarks, we have |]T0 



G7X = I d^xe'^^-'^ (vr^A(x)7r^A(0)) 

L (p' + Xi)(p' + X2)(p' + X3)K/3) 



+ XA ip' + XA?{V^ + M2J(J92 + M2)(p2 + M^, ) 

1 1 (p^ + Xi) (P^ + X2){P^ + Xs) 



(67) 



where A is a valence quark label, Xi normalized quark masses (with i = 1,2,3 for sea 
quarks), and M^^, M^, M^/ masses of neutral mesons in the sea sector (M^, = ml + 
Oixil^'o) )■ Clearly the limit mo — >■ oo exists, so the potential problem does not arise. 
The corresponding result in the quenched theory is 

'^AA — y^°) 

p^ + xa [p^ + XAr 

which shows that the mo — * oo limit cannot be taken, irrespective of the quark mass. 
The fact that one cannot integrate out $o in the quenched theory is consistent with the 
result, discussed in section ^ that $o should not be projected out of the theory, as it 
contributes to long range correlations. We observe again the central role played by the 
structure of str(TaTb) in the analysis. 



5.3 Reinterpreting previous calculations 

We have argued above that the $o field need not be included in PQ chiral perturbation 
theory because, in the massless limit, its propagator does not have a massless Goldstone 
pole, but instead has a pole at the rj' mass (if = 3). Since chiral perturbation theory 
does not converge for \p\ ^ IGeV, it will not, in general, be useful for \p\ ^ mrji. Thus the 
natural choice is leave $o out of the effective Lagrangian describing partially quenched 
simulations of QCD, perhaps adding it back as a mathematical device to simplify calcu- 
lations, as described previously in this section. 

Previous calculations in PQ chiral perturbation theory, e.g. Refs. |iy, |^, ^, |^, |23 



have, however, included $o as a physical field, rather than as a mathematical device. 
The reasons for this are partly historical (PQ chiral perturbation theory is an extension 
of quenched chiral perturbation theory, in which $o must be kept), and partly motivated 
by physics (there are regimes, e.g. large Nc or possibly = 2, in which $o is lighter 
than Aj^, and should be treated differently from other non-Goldstone hadrons — see also 
the discussion in Ref. 0). Our purpose here is to show that these previous calculations 
can be used to obtain the results in the theory without $0- 

We recall that including $0 as a physical field introduces several problems. First, 
chiral power counting is lost — adding a loop does not lead to an extra small factor of size 
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m^/A^, but rather to a factor of size mo/A^. If this latter factor is not small, diagrams 
with any number of $o loops must be included. Second, because of the anomaly, one can 
multiply terms in the effective Lagrangian by arbitrary functions of $o consistent with 
parity invariance. This introduces many new and poorly known parameters. Finally, 
the parameters ^i are those of unquenched chiral pertubation theory including the 77', 
which are related non-perturbatively, in a poorly known way, to those of the usual QCD 
chiral Lagrangian without the 77'. It is possible that one can mitigate some of these 
problems by assuming that the additional coupling constants are small (since they are 
1/iVc suppressed) .0 Such an analysis, however, inevitably involves assumptions which 
go beyond the systematic application of chiral perturbation theory. 

What we note here is that if one takes the limit mo — cxd in the results of these 
calculations, one recovers the results in the theory without $0 iii which the parameters 
are those of the QCD chiral Lagrangian without the r]' . The argument is a simple 
extension of that given in the previous section. If we expand in powers of $0) the total 
Lagrangian is 



of the previous section shows that they are irrelevant as mo 00. The remainder of the 
argument continues as before. 

What happens to the previous calculations when mo 00? First, the neutral propa- 
gator goes to an mo independent limit — the same limit discussed in the previous section — 
including, in general, light single and double poles. Second, "r^' loops" — those involving 
propagators with poles at m'^ ~ mo — give vanishing contributions. Third, all additional 
$0 couplings, such as a, do not contribute — as argued on general principles above. Thus 
the results of the calculations simplify considerably. 

In fact, this simplification has been noted and used previously, though not fully 
justified 0, |^, |^, It has been argued that, because mo ~ A^, terms proportional 
to m^/mQ are higher order in chiral perturbation theory and contribute only at the 
two-loop level. They can thus be dropped in NLO calculations. Contributions from 
T]' loops are not higher order (as noted above they can grow as In mo in dimensional 
regularization), but have been dropped by hand because they would not be present in 
QCD chiral perturbation theory calculations without the rj'. Together, this amounts to 
integrating out $0 perturbatively, whereas a non-perturbative treatment is required. Our 
results effectively provide such a treatment, and justify the ad-hoc procedures adopted 
previously. 



^•^Indeed, there has been significant success in applying quenched chiral perturbation theory to lattice 
data, despite the presence of the same problems (see, for example, Ref. p5[|). 



00 




(69) 
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6 Structure of the flavor- singlet propagator 



A signature prediction of PQ chiral perturbation theory at LO and NLO concerns the 
pecuhar features of the pole structure of correlation functions ("polology" for short). The 
effective theory predicts that propagators of flavor diagonal ( "neutral" ) mesons have both 
single and double pole singularities. In this section we show that, with some plausible 
assumptions about "normal" aspects of PQ polology, the occurrence of the unphysical 
double poles is a consequence of the symmetry structure of the theory. It follows that 
the double poles really are a feature of PQ QCD, and the grounds for this important 
prediction of PQ chiral perturbation theory are better established and understood. We 
first present the argument in terms of quark fields, so that it applies directly to PQ QCD, 
and then point out how the argument extends to the effective theory. 



6.1 PQ QCD 

For simplicity, we discuss only the case of degenerate sea quarks. 
Consider the two-point correlator of pseudoscalar quark bilinears, 

Gijkiix) = {Uij{x)llki{0)) , % = tr {QiQj-f^) . (70) 

The trace is over Dirac and color indices, and the order of Q and Q in the bilinears 
is chosen to simplify the transformation properties. U i = j and k = I, we call the 
correlators "neutral" , while if i ^ j and k ^ I we call them charged.]^ 

We first focus on the neutral propagator in a particularly simple theory, that with a 
single valence quark. A, its ghost. A, and one sea quark, S. In [Appendix E| we derive the 
constraints on G that follow from the graded vector symmetry together with Euclidian 
translation and rotation invariance. We find that the neutral propagator Gjjkk takes the 

following form in the {AA, SS, AA) basis: 



G = \ t u t \ . (71) 




The form holds at all separations and thus also in momentum space. By following quark 
propagators, we can interpret the elements of G as follows. (We use the language of 
perturbation theory, although the results hold non-perturbatively.) The off-diagonal 
term G^^^^ receives contributions only from disconnected diagrams, because the quark 
lines carrying the A flavor cannot be contracted with the A flavors. Gaaaa-, on the 
other hand, can get contributions from both connected and disconnected diagrams. All 
the disconnected graphs contributing to either of these terms in G have exactly the 
same structure, except for the interchange of A propagators with A ones. Since these 
propagators are equal, and there are no relative signs from Wick contractions, it follows 
that r in eq. (^) is the sum of all disconnected diagrams, and s is the sum of all 

"'^'*In the latter case, the unbroken vector symmetries require that i = I and j = k. This can also be 
seen by evaluating the correlator in terms of quark propagators. 
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connected ones.0 Similarly, t is the sum of all disconnected diagrams but with the quark 
line coupled to one of the external operators replaced by a sea quark (so that t = r if 
ruA = ms). Finally, u is the complete SS propagator, including both connected and 
disconnected contributions. Note that the relative signs of the various contributions can 
be determined simply by considering Wick contractions. 

Turning now to the more interesting case of > 1 degenerate sea quarks, it is 
straightforward to show that the same form of G holds for the neutral propagator in 
the basis in which SS is replaced by r]' cx "^j^i ^ SjSj. The quantities u, r, s, and t 
are of course different, with u being the t]' propagator. In this theory, there are N — 1 

additional neutral bilinears, but their correlators with AA, AA and rj' vanish since the 
latter are invariants under any symmetry transformation that involves only sea quarks, 
while none of the former is. Thus G is block-diagonal and we can consistently consider 
only the 3x3 block eq. (|7T|). 



For future use, we separate out from G the quantity 

Gc = diag{s,u,-s) . (72) 

We do this because we expect these quantities to have standard pole structure. The rj' 
propagator u is physical, and so has only single poles, with none of them being light. As 
for the connected propagator s, our key assumption in this section is that it has a light 
single-pole. This is what is predicted by PQ chiral perturbation theory at LO and NLO 
PI, and also what is observed in numerical simulations of PQ theories. 

In order to better understand the properties of the elements of G, it is useful to study 
the form of its inverse: 



G-^ = \ T U -T (73) 




with 



R=-^^-, S = i; r=-±; C/=i (74) 



us^ s us u 

UR-T^ 1 T 1 



t = — — ; u = —. (75) 



We now decompose G ^ as 



where 



US"^ ' S' US' U 



G-^ = + S , (76) 




(77) 



^^One can also show, using the graded symmetries, that s is related to correlators which are clearly 
charged: s = G^^^^ and s = Gabba if ^ha = tub (see Appendix E| ). 
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The form of eq. (|76D is the same as that in perturbation theory, with Gc corresponding to 
the free propagator, and E to the one-particle irreducible self-energy. Thus we interpret 
E as the self-energy contributions arising from disconnected diagrams involving at least 
some valence quarks or ghosts (disconnected contributions involving sea quarks having 
already been included). 

We can now state our second assumption: that the elements of S {R and T) are 
regular in momentum space at the position of the light pole in s. Our reasoning here 
is that S is one-particle irreducible with respect to the connected propagator, and thus 
does not have simple light poles. The only source for non-analyticity at the position of 
the light pole would be a deeply bound state of pseudo-Goldstone bosons, but this should 
not be present because the interactions are weak. 

Given our assumptions, we can now read off the pole stucture of G, by expressing it 
in terms of quantities with known, or assumed, behavior (simple light poles in s, around 
which u, R, and T are regular) :Q 

fs ~ Rs^ + T^us^ -Tus T^us^-Rs'^ \ 
G= -Tus u -Tus . (79) 

\ TV-Rs^ -Tus -s- Rs^ + TV/ 

Thus we see that there are light double poles (the terms), but no higher order poles. A 
more detailed analysis shows that the residues of the double poles vanish when = fns- 
This completes the argument for a single valence quark. What if there are two or more 
valence quarks? At the quark level, the addition of extra valence quarks does not change 
any of the elements of the block of G we consider. Thus the final form of G, eq. ([79D, 
is unchanged, and the same double poles are present. The argument holds separately 
for each valence quark, and thus Gbbbb has double poles of the same form as Gaaaa^ 
etc., although R and s depend on the mass of the valence quark. The argument can be 
generalized to study the structure of Gaabb (which has a double pole when = uib), 
but we do not give details here. 

6.2 PQ Chiral Effective Theory 

The previous argument concerned two-point functions of quark bilinears Iljj. In the 
derivation of the structure of G (eq. (|7ID), however, the Lagrangian of PQ QCD played 
no role, and the only properties of Iljj that were used were the transformation properties 
under vector symmetries. These properties are shared by the meson fields of PQ chiral 
perturbation theory, and therefore the propagators of these fields must also have the 
structure eq. (|7lD. Moreover, the statements about the analytic structure were based on 
the interpretation of the different components of G (the distinction between "connected" 

^^We note that can be simply inverted because the usual geometric series truncates after three 
terms: 

(Gc' + ^Gc- Gc^Gc + Gc^Gc^Gc ■ (78) 
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and "disconnected", "charged" and "neutral") which in turn used the language of Feyn- 
man diagrams and the tracing of quark lines. As discussed in section the effective 
theory can be formulated with $0? in which case the use of this language is still justified. 

It follows that in PQ chiral perturbation theory, with $o included, if the propagators 
for the charged mesons have only single poles, and the self energy function is analytic at 
low momentum, then there are double poles in the neutral propagators, and no higher 
singularities. Since the limit mo — > oo is well defined for these theories, and the low 
energy analytic structure is independent of it, we conclude that the above discussion also 
applies when the $o field is removed from the effective theory. 

Some of the assumptions and their implications are tested by the NLO calculation 
described in 0. There we found 

'7 1 



p2 + Ml^ ' ^ (p2 + Mls)/Zs + 5ss + ml ' 
R = ml + 5aa, T = ml + 6as- (81) 

Za and Zs are wave function renormalization factors and Maa and Mss meson masses, 
the expressions for which can be found in ||^. Relevant here is the fact that they are all 
independent of momentum and of mo, and that Mss and Maa are light. 6ss, ^aa, ^as 
are defined through 



Jab 



^LjXaXb + 48^(/ -Xa- Xb)\{Xa + Xb) ln(|(Xa + Xb))- (82) 



We see that the connected contribution to the propagator for the valence quark (s) 
has only a simple light pole, while the pole of the full sea quark 77' propagator {u) is 
heavy (at = — mg + ...). Also, the self energy components (i?, T) are clearly analytic 
in momentum. Finally, we collect all the terms into G, and take the limit mg — 00 to 

get 

-s I , (83) 




with 



r)2 _|_ Jij2 

e = -5aa + 25as - Sss - ^ \ (84) 



In summary, the argument presented above serves both as a confirmation that (given 
the validity of the assumptions) the pole structure of two-point functions predicted by 
LO and NLO PQ chiral perturbation theory is indeed that of PQ QCD, and an extension 
of that prediction to all orders in chiral perturbation theory. 

7 Conclusions 

The role of $0 in quenched and PQ chiral perturbation theory has been the main focus 
of this paper. We have shown that in order to reproduce the low momentum behavior 
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of two-point correlation functions of quenched QCD, $o niust be kept in the theory. On 
the other hand, in PQ QCD it does not give rise to long range correlations, in closer 
analogy to the r( in QCD, and should not be included. This point is key in carrying out 
the program outlined in [0, 0, |^, of using PQ simulations together with PQ chiral 



perturbation theory to determine the unknown constants that govern the low energy 
behavior of real QCD.|^ The central fact used in this program is that the parameters 
of the chiral Lagrangian in QCD (with 2 or 3 light flavors) and in PQ QCD (with the 
corresponding number of sea quarks) are the same. In the presence of however, the 
PQ chiral Lagrangian matches the unquenched chiral Lagrangian in which the 77' field is 
present. The latter theory does not have a low energy expansion (for the physical values 

= 2,3 and A^coior = 3), and its relation to low energy hadron phenomenology cannot 
be calculated perturbatively. What we have shown is that this problem does not arise 
because the PQ chiral effective theory can and should be formulated without $o- 

In seeming contradiction to what has just been stated, there are technical benefits 
from keeping $0 iii the effective theory. In this paper we have shown how $0 can be 
included as an auxiliary field with a mass term ^0*^05 ^^^1 its effects can be then removed 
by taking mg to infinity. This also establishes the status of previous results in PQ chiral 
perturbation theory in which $0 was kept. By taking mo 00, all effects of $0 are 
removed from these results, irrespective of whether other $0 couplings were included. 

The role of $0 is tied in with the more general theme of the foundation and justification 
of PQ chiral perturbation theory. We have addressed this issue by attempting to repeat 
the line of reasoning that leads to the standard chiral Lagrangian. As a first step, we have 
identified the full symmetries of PQ QCD. We then argued that the symmetry breaking 
pattern in this theory can be derived from the symmetry breaking pattern of QCD. 
Goldstone's theorem, with the use of the appropriately generalized Ward identities, then 
leads to the conclusion that two-point correlation functions of operators associated with 
generators of broken symmetries have low-lying poles. We discussed in some length the 
construction of the effective theory for the fields that have the same quantum numbers 
as these operators. This theory is guaranteed to recover the low energy behavior of two- 
point functions of the chiral currents and densities. In the absence of a Hilbert space, 
however, we do not know how to show that the long range behavior of general n-point 
functions can be attributed to the singularities of the two-point functions. This is a 
crucial implicit assumption that is made when one uses the PQ chiral effective theory. 

Lacking a general argument to justify PQ chiral effective theory, we have focused 
instead on one of the strikingly unphysical aspects of PQ chiral perturbation theory, the 
existence of light double poles in propagators of flavor-neutral mesons. We have demon- 
strated that the existence of these double poles (and the absence of higher singularities) 
follows from the assumption that the propagators of charged mesons have only simple 
poles.Q The proof involves only the symmetries of the theory, symmetries that are shared 
by the underlying microscopic theory and the low energy effective theory. We learn two 
things from this result. First, that this unphysical feature of the effective theory is cor- 

^^First results from this program have recently been presented pTf . 

-"^^This assumption can itself be derived in the massless theory (as in section ||), but not in the inter- 
esting case of massive quarks with mv 7^ 
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rectly representing the properties of the underlying PQ theory. And, second, that the 
pole structure seen in LO and NLO chiral perturbation theory will hold also at higher 
orders — there will only be single and double poles. 

Finally, we note that an interesting consistency check of our results can be obtained by 
taking the valence quark masses to be much smaller than the sea quark masses, though not 
so small that enhanced chiral logarithms [§, proportional to mslnmv, invalidate chiral 
perturbation theory.^ In this regime the PQ theory has a "light" sector, with correlators 
having poles at M^^^^ oc niv, and a "heavy" sector with poles at M^^^^^ oc mg. We expect 
the relevant degrees of freedom in the light sector to be the valence quarks and ghosts 
alone, and thus that it should be described by an effective quenched SU{Nv\Nv) chiral 
Lagrangian. In particular, this Lagrangian should contain the quenched $^ field, despite 
the absence of the $o in the underlying PQ chiral Lagrangian. We also expect that 
additional terms, such as a should appear in the quenched effective Lagrangian. 

This issue can be investigated analytically, since both valence and sea quarks are in the 
chiral regime. We have checked that our expectations are indeed borne out, by matching 
pole positions at one-loop obtained from the underlying PQ theory and the effective 
quenched theory. We find, for example, that a PQ theory having degenerate sea 
quarks matches with a quenched theory having ttIq = 2Boms = Xs a = 1 (as well as 
small, ms and dependent, shifts in Bq and / between the two theories). 
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Appendix A Real and Fake Symmetries 



In this appendix we discuss the symmetry group of the PQ theory defined in eq. ([l|), and 
the resulting Ward identities. We show that, ahhough the fiavor symmetry group differs 
from the naive expectation SU{Nv + N\Nv)l(^ SU{Nv + N\Nv)r, the Ward identities 
coincide with those derived assuming this "fake" symmetry to hold. This appendix is 
based in part on the analogous development for the quenched theory worked out in 



Ref. 17 



Quark Sector Symmetries 

Consider first only the quark part of the action, eq. 



{^L^^pL + i'RlP^R + ^im^R + ^Rm^L) (85) 
In the massless limit it is invariant under transformations of the form 

i^L,R GL,Ri>L,R, i^L,R ^ tpL,RGl% (86) 

where Gl and Gr need only be non-singular, Gl r E GL{Ny + A^). Requiring that the 
functional measure be invariant reduces the symmetry to 

SL{Ny + N)l ® SL{Nv + N)r ® U{1) . (87) 

In the following we focus on the fiavor symmetries, and do not show the overall U{1) 
phase symmetry. 

The group ( 187|) is larger than the symmetry group of the Hamiltonian, SU{Nv + 
N)l <S) SU{Nv + N)r, because ip ^ 'ip'^'Jo in the functional integral formulation. To 
understand the physical significance of the enlarged symmetry, we consider the result- 
ing Ward identities. The infinitesimal transformations take the form (choosing Gl for 
illustration): 

Gl = exp(aT) ^ 1 + aT (88) 
6tpL = (yTi)L, S^pL = -ai^LT, (89) 

where T is an arbitrary traceless, hermitian {Ny + N)"^ matrix, and the small parameter, 
a, is complex. The usual, unitary, symmetry transformations correspond to a being 
pure imaginary, and thus have half as many parameters. Choosing a to be space-time 
dependent, we have (still in the massless theory) 

J JlJ = i>Ll,Ti^L, (90) 

0{y)^0{y) + a{y)5'^P0{y), (91) 



29 



where (9 is a local operator.^ Such a transformation can be seen as a change of variables 
in the functional integral, with unit Jacobian, and therefore leads (in the example of a 
single operator) to 

= 6 {0{y)) = J a{x)d, {j^^Ji^My)) dx + a{y) (d^Oiy)) . (94) 

Note that eq. (|9^) contains a but not a*. Thus, although it appears that two equations 
can be obtained for each independent matrix T (for real and imaginary a), in fact only 
one Ward identity results: 

(jl?(^)0(y)) = -5ix - y) {5^P0{y)) . (95) 

Thus the unitary subgroup (imaginary a) is sufficient to generate all the Ward identities 
implied by the full symmetry group. 

When the masses are not zero, another term is added to eq. (|5^), corresponding to 
the variation in the action because of the mass terms. These terms, however, like the 
operator O, do not involve complex conjugates of the quark fields, and therefore their 
variation under the transformation involves only a. As in the massless case, this implies 
that there is only one Ward identity (a modification of eq. (^)) corresponding to each 
generator T. 

Clearly the same analysis holds for the right-handed transformations. We conclude 
that the unitary subgroup, SU{Nv + A^)l ® SU{Nv + N)r, is sufficient to generate all 
the Ward identities implied by the full symmetry group. 



Ghost Sector Symmetries 

In the PQ QCD partition function, eq. (0), one integrates independently over the Grass- 
mann fields and ip. The integral over the commuting ghost variables, however, con- 
verges only if it has a Gaussian structure: 

J DxDx'' exp{—x'' Ax) , (96) 

where the hermitian part of A must be positive definite. Since Ip is anti-hermitian, this 
constraint applies to the mass term. Thus in the ghost part of the action, 

j {i^JPi'L + i^RlPi^R + ^L^i^R + ^Rm^JL) , (97) 

^"^It is convenient to use a somewhat inconsistent notation in which the variations of fundamental 
fields and generic operators are defined differently. Thus, in the case of left-handed transformations of 
the quark fields 

+ 5l^, (92) 
0^0 + uSlO, (93) 

so that, in fact, 6lO is not an infinitesimal quantity. 
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we must identify ■^^ = ■^jj and ijjji — ipy. The ghost part of Sp is then: 

j ^\lp^R + ^^Ip^L + ^{m^L + ^^m^R ■ (98) 
The symmetries of the kinetic terms alone are thus 

^L^G^L, ^r^{G-')Hr, (99) 

where G G GL{Nv). The anomaly reduces this GL[Nv) symmetry group to a product 
of SL{Nv) and an overall phase rotation. 

Ward identities are derived from infinitesimal local transformations, G — ex.p{aT) ~ 
1 + aT, with T a traceless, hermitian, Ny x Ny matrix, and a complex, leading to 

SipL = aTipi, , 6ipR = -a*T2pR , 6ipl = a*iplT , 6ip\^ = -aip\{r . (100) 
Local operators, 0{7jjL,''pL,'4^R,i'R), transform like 

SOiy) = (^T^PL - i^^RT^] - ( ^Ti.R - ^'lT^] ■ (101) 

" V ' " V ' 

=5lO =SiiO 

One obtains (for the case of the expectation value of a single operator) 

^ a{x) {djPix)Oiy)) - a{xr {dg\x)0{y)) = -a{y) (^^©(y)) + a(y)* {5RO{y)) , 

(102) 

with 

& = AL^,Ti^L,R. (103) 

By taking a real and imaginary the full set of Ward identities are seen to be equivalent 
to the equations 

(dj^Ri^My)) = -S{x - y) {Sr.,RO{y)) . (104) 

The generalization to m 7^ is straightforward. 

The resulting identities are exactly those that one would obtain were one to pretend 
that ip and ip were independent, and that the symmetry group was SU (iVy) l^SU {Nv)r. 
We note, however, that, unlike the situation in the quark sector, the true symmetry group 
SL{Nv) does not contain the "fake" symmetry group SU{Nv)l® SU{Nv)r. While they 
do share a common vector SU (Ny) subgroup, and have the same number of generators, 
the axial transformations (left and right handed fields rotating oppositely) differ. We 
also note that the use of complex a is crucial for obtaining all the independent Ward 
identities. 
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Graded Symmetries 

To complete the symmetries we need to consider graded transformations which rotate 
Grassmann and ghost fields into each other. Once we do this, we are necessarily con- 
sidering the ghost fields to be commuting elements of a Grassmann algebra, with (in 
the notation of Ref. [|l^]) both a "body" or "base" — the usual complex scalar field — 
and a "soul" — composed of products of an even number of Grassmann fields, and thus 
nilpotent. The following question thus arises: What constraint does the requirement of a 
convergent functional integral impose on ip and ipl We argue at the end of this Appendix 
that the answer is 

I body = I body ' I body = \ body ' 

i.e. the relations discussed above hold for the bodies of these quantities, but not for the 
souls. In this subsection we discuss the consequences of this constraint on the symmetries. 
To do this, we return to the notation 

QIr = ^Ib) > Ql,R = {^L,R, k,R) ■ (106) 

The most general anomaly-free fiavor transformation which leaves the kinetic term 

J Ql^Ql + QrIPQr, (107) 

invariant is 

Q_L _LQl Q_r _RQr ^-j^gg^ 

Here L and R are independent SL{Ny + N\Ny) matrices, except that they must maintain 
the constraint (|105| ). If we write 



Lqq Lqg 



(109) 



(and similarly for L ^, R and R ^) to denote the quark-quark block, quark-ghost block, 
etc., of the different matrices, then the constraint is 

^sslbody = i^gg) Ibody- (110) 

In deriving this, we have used 



(-^ ) 99 1 body (-^99 1 body) "^99 I body (m) 

and related results, which follow from the fact that only the product of "bodies" con- 



tribute to the body of a product. The matrices L and R satisfying the constraint eq. ( |110| ) 
form a subgroup of SL{Nv + N\Nv) l ® SL{Nv + N\Nv)r. The exponential parameter- 
ization of elements of this subgroup is 

L = exp(z$i), i? = exp(i$R), str (<l>i,^) = , (^l)^^ | [^.jy = (*«)99 1 body ' (l^^) 
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To derive Ward identities we consider infinitesimal transformations of the form in 
eq. (|112|) . Transformations in the quark-quark and ghost-ghost blocks lead only to the 



identities described in the previous subsections. In particular, while the full symmetry 
group has axial transformations in the ghost-ghost block which were not considered above 
(in which ($L,fi)g3 are both pure soul) these do not lead to independent Ward identities. 

Additional Ward identities do arise from purely graded transformations. These are 
derived by considering i^l.r = oil,rT-, with T hermitian matrices contained entirely in 
the quark-ghost and ghost-quark blocks, and ai^R anticommuting parameters. Note that 



there are no constraints on aL,R from eq. (|112D , and so the derivation of Ward identities 



follows the same steps as in the quark sector, except that one must keep track of the 
anticommuting nature of a. The result is a single independent identity of the form of 
eq. (p5| ) for each "off-diagonal" generator T. As in the quark sector, the identities are 
the same as those that follow from the unitary subgroup in which ^l,r are constrained 
to be hermitian. The extra freedom of complex parameters does not lead to additional 
identities. 

Combining the results from all infinitesimal transformations, we see that all the Ward 
identities could have been obtained if one had assumed the fake symmetry group SU{Ny + 
N\Nv)l ® SU{Nv + N\Nv)r. 

Vector and axial transformations and currents 



The VEV of eq. ( |T5D (and also the mass term with M oc /) breaks the chiral symmetry 
of eq. ( |108D down to its vector subgroup: 

L = R^Ve SL{Ny + N\Ny) , V,,l^^^ = {V-^') . (113) 

The corresponding "fake" group, sufficient for deriving vector Ward identities, is the 
subgroup V G SU{Nv + A^|A^y), for which the constraint in eq. ( p.l3| ) is automatically 
satisfied. 

The axial transformations, the generators of which are broken by the VEV, are given 

by 

L = R-'^Ae SL{Ny + N\N), A..^^^^ = (A,,) . (114) 

Here the fake transformations have A G SU{Ny + A^|A^), and are not contained in the 
transformations of eq. ( |114| ) because the constraint is not satisfied. 



By combining infinitesimal left- and right-handed transformations in the appropriate 
way, we can derive vector and axial Ward identities. They take the same form as eq. (pSj), 
with L — > V, A, and contain 

Jv,U^) = QL-f,TQL{x) ± QRl,TQn{x) (115) 

and 

€l0 = 6r0±6'^P0. (116) 
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The only subtlety here is that, when deriving Ward identities for graded transformations, 
the factor of a{y) which is pulled out of eqs. (|90| - p4D should be accompanied by the sign 



which results when moving a past Q. This is needed to be consistent with the definition 
of and impacts the definition of 5^"^^. 

Convergence considerations 

First recall that integrals over real c-numbers lead to results of the same form as over 
ordinary numbers [|18| 

f dqf{q)=F{h)-F{a), F' = f . (117) 

J a 

Here functions of c-numbers are defined by Taylor expansions, 

oo 

m = f {iL.y) + E (^Loui)" /^"^ (^Ibod.) ' (118) 

n=l 

where the sum actually truncates because q\^^^^ is nilpotent. In our case we are interested 
in Gaussian integrals: 

/(g) = exp(-mg2) , a = -oo + aj^^^j , b = +oo + b\^^^^, (119) 
with Re(m) positive. Since /*■"■' (±oo) = for all n, it follows that eq. (117) is actually 



independent of the souls of a and b. Similarly, if we change variables by a quantity which 
is pure soul, q' = q + 5q, Sq\^^^^ = 0, we do not have to change the limits of integration.]^ 
This is true for an arbitrary change in soul — in particular, it does not need to be real. 
Now consider a two-dimensional Gaussian integral 

/ + 00 
dqidq2 exp [-m{ql + qj)] . (120) 
-oo 

As we have just seen, the bodies of qi^2 are real, but their souls can be arbitrary without 
changing the value of the integral. Thus if we change to "complex" variables 

q = q^ + iq2, q = q^ ~ iq^ , dqdq = dqidqz , (121) 

then the integral takes on the usual complex Gaussian form 

dqdq exp {—mqq) , (122) 



except that q = q* holds only for the bodies, and not for the souls (since the souls of gi_2 
are not real). The integral itself has the value given by eq. ( |12CI| ), i.e. ir/m. Note that 
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Note that this is not true for general integrals: one must take care when making changes of variables 



involving nilpotent parts, since they can lead to so-called anomalies 11, 19 
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m can also have an arbitrary soul — what matters for convergence is that the real part of 
its body is positive. 

The generalization to many complex variables is straightforward. Let q and q now 
represent a vector and transposed vector, respectively. Then the integral 

J dqdqexp{-qMq) (123) 

is convergent if M (taken to have no soul for now) is hermitian and positive, and if q = q^ 
holds for the bodies. To see this we diagonalize M: M = U'^DU, with D being diagonal 
and positive, and U unitary. Thus if we use the variables q' = Uq and q = qW (which 
leaves the measure unchanged), then the integral factorizes into product of integrals of 
the form of eq. ( [I22| ), each of which is convergent as long as ^'jj^^^y = "y'lbo^y' '^^^^ 
relation is maintained by the unitary transformation back to the unprimed fields. The 
argument goes through if M has an arbitrary soul, since the integrand can be expanded 
in powers of this soul, and each term is convergent. Note that the general symmetry 
transformation, eq. ( |1U8| ), maintains the hermiticity of the body of M. 



Appendix B The diagonal generators of U {Nv+N\Nv] 



Here we collect some useful results concerning the generators of graded symmetry groups. 
These are represented by the hermitian {2Nv + A^)^ matrices labeled T in the foregoing. 
Note that the same generators serve for both the true and fake symmetry groups. We 
consider here the properties of the diagonal generators. 

Let Aa be the Ny + — 1 diagonal generators of SU{Nv + N), chosen to satisfy 
tr(AaAft) = 6ab- Similarly, let Aa be the A^'y — 1 diagonal generators of SU{Nv), normalized 



in the same fashion. We define T„ 



N + 2A'"v' — 2 to be the set of matrices: 



Ny + N quark indices{ 
Ny ghost indices! 



where we use a schematic notation. They satisfy 

str(r,) = 

and 



(124) 



tr(Aa) 
-tr(Aa) 



(125) 



str(raTfe 




(126) 



Two more generators need to be defined to complete the basis for the diagonal part 
oiUiNv + N\Nv). 
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Partially quenched (A^ 7^ 0) 

We choose 



2Nv+N-l 



2^yNNv{Nv+N) 



{NI - {2Nv + N)I) 



and 



where 



T2NV+N - -^I, 







(127) 



:i28) 



(129) 



The last element, T2Nv+n, has non-vanishing supertrace, and generates the anomalous 
f/(l) factor of U{Nv + A^|A^y). 

Considering now all 2Nv + N generators, it is straightforward to check that 



{str{T,n)) = {gab) = diag(l. 



,1,-1, 



-1,-1,1). 



(130) 



Nv+N-l Nv-l 



Quenched (A^ = 0) 

In this case we choose 



T2NV-1 — ^ 
T2NV = JFh^^- 



(131) 
(132) 



While the identity is straceless, / is not, and is taken as the generator of the anomalous 
f/(l) factor. Eq. (|130|) becomes 



/I 



{sti (TaTb)) = (gab) 



\ 



V 



1 

1 0/ 



(133) 



Appendix C Another argument concerning (^q 

In this appendix we give an alternative argument concerning the status of $0 in quenched 
and PQ chiral perturbation theory. We consider two-point functions of pseudoscalar 
densities. 



G'(^'^')(x) = (0(^)(a;)0(^')(O)) 



(134) 
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where T, T' run over all the generators of U{Nv + N\Nv) (and thus include the identity). 
For the diagonal generators, we use the basis given in [Appendix B[ We assume that all 



quark and ghost masses are equal, although they do not have to vanish.^ 

Consider first the PQ theory. As long as > 2, some of the generators lie entirely 
in the sea-quark sector. Choosing T and T' of this form, we can use chiral perturbation 
theory for QCD-like theories to infer that the resulting correlator has a pseudo-Goldstone 
boson (PGB) pole if T = T', while flavor symmetry implies that it vanishes if T 7^ T'. 
Thus we know that 

G^^'^'\x) = 6t,t'Gpgb{x) , T, T' e T,ea . (135) 

We note for future reference that G'pGB(a;) is "connected" in the sense that the only 
contractions which contribute are those in which the two bilinears are connected by 
quark propagators. 

We can extend this result to all the straceless generators using the graded vector 
symmetry (which, as argued in the text, is not spontaneously broken), with the result 

G'-^'^'^x) = stT{Tr)GpGB{x) str(T) = str(r) = . (136) 



This can be shown either directly using the symmetry, along the lines of [Appendix E| , or 



by a direct comparison of the contributing contractions. In the latter case, the overall 
factor str(TT'), which can be of either sign (see Appendix B ), accounts for the signs 



arising from fermionic Wick contractions. The result ( [136[ ) shows that there are PGB 
poles in the correlation functions for each of the generators of SU{Nv + N\Nv). This 
agrees with the result obtained in and implies that the corresponding fields should 
be included in the effective Lagrangian. 

Now consider the correlation functions involving the interpolating field corresponding 
to $0; which are obtained by setting T and/or T' to the identity. Flavor symmetry implies 
that if T = J and str(T') = 0, or vica-versa, then the correlator vanishes. If T = T' = /, 
then, by inspecting contractions, one can show that the correlator is proportional to that 
for the 77' [p^ : 



G(^'^)(x) = ($o(x)$o(0)) oc {v'{x)r]'{0)) (137) 

where the r]' is the flavor singlet meson in the sea sector. Valence quark and ghost 
contributions completely cancel. Now, due to the axial anomaly, we know that the 77' 
correlator does not have a light pole. Thus none of the two-point functions involving 
$0 have PGB poles. Conversely, to describe the long distance parts of the two-point 
functions in the PQ theory we do not need to consider correlators involving $o- 

We now contrast this analysis with that for the quenched theory. By examining 
contractions, or using the graded flavor symmetries, one can show that 

Gg^'^'^(x) = str(TT')Gconn(a;) + str(T)str(T')Gdisc(2;) , (138) 



^^Note that the PQ theory differs from the unquenched theory (with N flavors) even if aU masses are 
equal. This is because the extra fields in the PQ theory allow one to separately determine certain Wick 
contractions which always arise in certain linear combinations in the unquenched theory. 
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where "conn" and "disc" refer to connected and disconnected contractions, and T and 
T' run over the generators of U{Nv\Nv)- Since there is no sea sector in the quenched 
theory, we cannot rely on experience with QCD to imply that there are PGB poles in some 
channels. Instead, based on numerical data, we assume that there is such a pole in the 



connected correlator Gconn{x)- Then we see from ( p.38| ) that, if we use the basis explained 
in [Appendix B| , there is a PGB pole in correlators corresponding to all the generators of 



SU{Nv\Ny) except T = I. This if because, for such generators, str(T^) 7^ (so the first 
term in eq. ( |138|) is present) but str(T) = (so the second term is absent). However, for 
T = T' = /, both str(T) and str (TT') vanish, and so Gq'^\x) = 0. On the other hand, 
if T = / and T' = I (the anomalous generator), or vica- versa, then str(TT') 7^ 0, and so 

G'g''^CxGconn(x). (139) 

Thus this cross-correlator also has a PGB pole. Finally, if T = T' = /, then 

Gg-'^ oc Gdisc(x) , (140) 

about which the present analysis says nothing, although quenched chiral perturbation 
theory predicts a double PGB pole. 

The important conclusion is that, to include all channels which have PGB poles in 
them one must include both (f)^^^ and (p^^^ oc $o- 



Appendix D G in the limit mo oo 

In this section we calculate the propagator G(o-) in the large rriQ limit, where the inverse 
propagator is given by 

with A, B, C, and d independent of rriQ. Note that A is a square matrix, B and C 
a column and row vectors, respectively, and d a number. Since d appears only in the 
combination ttLq + ci it can be dropped in the large rriQ limit. 
We write 

G,., = {c, . (142) 

To learn about the tuq dependence of the blocks of G(cr) we consider the equation 

GmG(-4=/. (143) 
One of the equations contained in eq. ( |143| ) is 

C'A + d'C = 0, (144) 
which implies d' ^ C, where the tilde refers to the scaling with mo. 
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Another equation in eq. (|143|) is 

C'B + d'ml = 1. (145) 

Since d' and C scale the same, only the second term on the left hand side of eq. ( |145| ) is 
important when mo — *■ oo, and we conclude that 

d',C'^^. (146) 



rriQ 



Similarly, from 



we get 



which is then used in 



A'B + B'ml = (147) 



B' ~ -\a', (148) 



m, 







A'A + B'C = 1 > A'A=1. (149) 

mo^oo 



In the last equation we see that cannot be inverted when nio oo unless A is a 
non-singular matrix. 

Putting everything together, in the large mo limit 

A' = A-' + o(^] =0(1), (150) 
B',C',d' = o(^), (151) 



mo 



and therefore 



mo 



Appendix E The Structure of the Propagator from 

Graded Symmetries 

In this appendix we derive the constraints on the structure of the pion propagator that 
imply the form eq. (ffll). The following symmetries are used: 

• Independent phase rotations of individual flavors. These form a subgroup of the 
vector transformations, eqs. (|108[ - |113| ), where one takes 

= diag(exp 6^^, exp 6^5, exp 6*^), (153) 

and 6a, ds and 6^ are independent. Under a phase rotation of only the flavor m, 

Gijkl Gijkl exp(6'm(5im — Sjm + ^km — ^Im))- (1^4) 
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SU{1\1) transformations that involve only A and A. These too are vector transfor- 
mations, this time with 



(155) 




where 

U 



a b 

c 



J) e SU(1|1). (156) 



• Gijki = {—)Gkiij- This symmetry follows from 

(n,, (x) liki (0))= (n,, (-x) Hh (0)) 

= (n,,(0)nH(a;)) (157) 
= (-) {Yiu (x) n,, (0)) , 

where the first equation is obtained by rotation, the second by translation, and the 
(— ) sign in the third is only needed when both Iljj and 11^; are fermionic fields. 

The invariance of G under transformations of the form eq. ( |154] ) implies that the in- 
dices must be paired up (quark lines must be followed, corresponding to legal contraction 
of quark fields). The non- vanishing elements of G therefore take the form Gujj or Gijji. 
The implications of the SU symmetries are slightly less straightforward. We first 
form the following 2-indexed objects out of the elements of G: 

G^nk{x) = ((n(x)n(o)),,) 

j=A,A 

Yl ^i^jjikix) = (str(n(x))nifc(0)) ^^gg^i 

j=A,A 

Gssikix) = (n5s(x)n,,(o)) 

Gisskix) = {U,six)Uskm ■ 

The matrix notation (11 with no indices, matrix multiplication, and the strace symbol) 
refers to 2 x 2 matrices in the A — A subspace. In a consistent manner, i,k ^ {A, A}. 
All of the combinations above transform similarly under SU{1\1) transformations: 

m,n=A,A 

or, using the 2x2 notation again, 

O UOUl (160) 
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Since SU{1\1) is a symmetry, each of the combinations are invariant under these transfor- 
mations. Direct examination shows that this imphes that each O must be proportional 
to the identity. From this argument we obtain the following set of equations: 



Gijjk = r6ik ^ < _ (161) 

^ { '-'AAAA "T '-^AAAA ~ ' 

E ^,G„., = ^ ( ^-"^^ : ^--^ : " (162) 

. . 7 I '-^aaAa '-'aaaa — ^ 

]=A,A V 

GsSik = tSik =^ GsSAA = Ggg^^ = t (163) 

Gissk = v5ik =^ Gassa = G^gg-^ = v. (164) 

Finally, with the use of eq. ( p.571) , these equations can be solved to yield: 

Gaaaa = r + s (165) 

Gaaaa = r-s (166) 

Gaaaa = ^AiAA = (167) 

GaaaA = -Gaaaa = ^ (168) 

GssAA = Gaass = Ggg^^ = G^^gg = t (169) 

Gassa = Gsaas = G^55A = -GsAAs = (1^0) 



The last independent element of G is Gssss- The form shown in eq. (^) follows (what 
appears there is the restriction of G to the subspace of AA, SS and AA). 
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